Journal  of  Power  Sources  196  (2011)  5685-5696 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


in 

SllbaidtS 


Reciprocating  air  flow  for  Li-ion  battery  thermal  management  to 
improve  temperature  uniformity 

Rajib  Mahamud,  Chanwoo  Park* 

Department  of  Mechanical  Engineering,  University  of  Nevada,  1664  N.  Virginia  Street,  Reno,  NV  89557-0312,  USA 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  13  January  2011 

Received  in  revised  form  14  February  2011 

Accepted  17  February  2011 

Available  online  26  February  2011 


Keywords: 

Li-ion  battery 
Thermal  management 
Reciprocating  flow 
Electrical  vehicle 
CFD 

Lumped-capacitance  thermal  model 


The  thermal  management  of  traction  battery  systems  for  electrical-drive  vehicles  directly  affects  vehicle 
dynamic  performance,  long-term  durability  and  cost  of  the  battery  systems.  In  this  paper,  a  new  battery 
thermal  management  method  using  a  reciprocating  air  flow  for  cylindrical  Li-ion  (LiMn204/C)  cells  was 
numerically  analyzed  using  (i)  a  two-dimensional  computational  fluid  dynamics  (CFD)  model  and  (ii)  a 
lumped-capacitance  thermal  model  for  battery  cells  and  a  flow  network  model.  The  battery  heat  gener¬ 
ation  was  approximated  by  uniform  volumetric  joule  and  reversible  (entropic)  losses.  The  results  of  the 
CFD  model  were  validated  with  the  experimental  results  of  in-line  tube-bank  systems  which  approx¬ 
imates  the  battery  cell  arrangement  considered  for  this  study.  The  numerical  results  showed  that  the 
reciprocating  flow  can  reduce  the  cell  temperature  difference  of  the  battery  system  by  about  4  °C  (72% 
reduction)  and  the  maximum  cell  temperature  by  1.5  °C  for  a  reciprocation  period  of  r=  120  s  as  com¬ 
pared  with  the  uni-directional  flow  case  (r  =  oo).  Such  temperature  improvement  attributes  to  the  heat 
redistribution  and  disturbance  of  the  boundary  layers  on  the  formed  on  the  cells  due  to  the  periodic  flow 
reversal. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electric  vehicles  (EVs)  and  hybrid  electric  vehicles  (HEVs)  are 
considered  to  be  sustainable  and  environmentally  friendly  trans¬ 
portation  options  in  the  midst  of  rising  oil  prices  and  global 
climate  changes  due  to  greenhouse  gas  emission.  The  traction 
battery  systems  for  the  electrical-drive  vehicles  houses  a  large 
number  (hundreds  to  thousands)  of  high  energy-density  recharge¬ 
able  batteries  such  as  Li-ion,  and  Ni-MH  to  satisfy  the  demanding 
vehicle  power  requirements.  The  battery  system  performance 
greatly  affects  all  the  vehicle  dynamic  and  electric  parameters 
such  as  acceleration,  regenerative  braking  and  battery  charg¬ 
ing/discharging. 

The  battery  thermal  management  (including  both  cooling  and 
heating)  is  essential  for  electrical-drive  vehicles  to  maintain  proper 
cell  temperature  conditions  in  all  weather  conditions.  The  thermal 
imbalance  among  the  cells  significantly  affects  short  and  long¬ 
term  performances  of  the  vehicle  battery  systems.  According  to 
the  Arrhenius  law  of  battery  electrochemistry,  the  battery  reac¬ 
tion  increases  exponentially  with  the  battery  cell  temperature.  As 
a  result,  hotter  cells  degrade  more  quickly  than  colder  cells  [1-4] 
and  those  few  overheated  cells  result  in  shortening  the  lifetime  of 
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a  whole  battery  pack.  The  lifespan  of  the  Li-ion  cell,  for  example, 
is  reduced  by  about  two  months  for  every  degree  of  temperature 
rise  in  an  operating  temperature  range  of  30-40  °C  [5].  The  Li-ion 
battery  system  is  required  to  maintain  the  maximum  cell  temper¬ 
ature  below  40  °C  and  the  cell  temperature  difference  below  5°C 
for  a  full  lifespan  [2]. 

The  primary  objectives  of  the  battery  thermal  management 
are  (1)  to  limit  the  battery  cell  temperatures  below  the  allowed 
maximum  cell  temperature,  Ts,max;  (2)  to  minimize  the  cell  tem¬ 
perature  difference,  A Ts;  and  (3)  to  maintain  the  cell  temperatures 
within  the  operating  range  for  optimum  performance  and  longevity 
of  the  battery  system.  Air  (ambient  and/or  air-conditioned)  has 
been  the  coolant  of  choice  because  of  reliable  and  simple  bat¬ 
tery  cooling  systems.  Specially,  the  battery  cooling  would  be  very 
challenging  in  the  summer  (ambient  air  temperature  of  50  °C) 
without  using  vehicle  air-conditioning  system.  The  conventional 
battery  cooling  systems  utilize  the  one-directional  air  flow  mov¬ 
ing  from  the  inlet  and  outlet  of  the  battery  cooling  systems.  As 
a  result  of  the  convective  heat  transfer  along  the  air  stream,  the 
air  temperature  at  the  outlet  is  always  higher  than  the  inlet  and 
the  downstream  cells  near  to  the  outlet  are  likely  hotter  than  the 
upstream  cells.  Most  of  the  conventional  battery  thermal  manage¬ 
ment  systems  actually  use  flow  distributors  (baffle),  which  also 
serve  as  battery  holder,  to  create  a  strong  convective  heat  transfer 
condition  for  the  downstream  cells  in  an  attempt  to  compensate 
the  unfavorable  air  temperature  for  the  convective  heat  transfer 
[3]. 
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Nomenclature 

A  area  [mm2] 

Bi  Biot  number 

C  correction  factor 

cp  specific  heat  capacity  [J  kg-1  K-1  ] 

D  diameter  of  a  battery  cell  [mm] 

dEoddt  entropic  coefficient  [mV K-1  ] 

E  cell  potential  [V] 

/  friction  factor 

h  convective  heat  transfer  coefficient  [W  m-2 I<-1  ] 

/  electrical  current  [A] 

k  thermal  conductivity  [W  m-1  K-1  ] 

L  length  of  a  battery  cell  or  characteristic  length  [mm] 

l  duct  length  [mm] 

Mf  mass  flow  rate  [kgs-1] 

n  total  number  of  tubes  (battery  cells)  or  index  num¬ 

ber  of  the  last  tube  (or  battery  cell) 

NTU  number  of  transfer  unit 

Nud  Nusselt  number  based  on  the  cell  diameter  D 

NuD  c  corrected  Nusselt  number 

Pr  Prandtl  number 

p  pressure  [Pa] 

Q  heat  transfer  rate  [W] 

R  thermal  resistance  [K  W-1  ] 

Re  Reynolds  number 

Re  internal  electrical  resistance  of  a  battery  cell  [m£2] 

5  spacing  between  cells  [mm] 

Sr  battery  heat  generation  [W] 

T  temperature  [I<] 

T  time-averaged  temperature  [K],  f  =  f  Tdt/  f  dt 

t  time  [s] 

U  velocity  [ms-1] 

V  volume  [m3] 

V  volumetric  flow  rate  [m3  s-1  ] 

w  width  of  the  computational  domain  [mm] 

Greek  _ 

AT  time-averaged  temperature  difference  [K],  AT  = 

Ts,  max  -  Ts  m[n 

s  porosity 

fi  viscosity  [Pas] 

p  density  [kg m-3] 

r  reciprocation  period  [s] 

Subscripts 
b  battery 

c  cyclic  or  characteristic  or  correction 

D  diameter 

e  electrical 

/  fluid  phase 

i  battery  node  index  or  duct  inlet 

irr  irreversible 

ku  surface  convection 

L  longitudinal  direction  or  duct  outlet 

LMTD  log  mean  temperature  difference 

max  maximum 

min  minimum 

n  index  number  for  the  last  battery  cell 

oc  open  circuit 

o  duct  outlet 

p  pressure  or  particle 

r  battery  reaction 


rev  reversible 

s  solid  phase  or  surface 

T  transverse  direction 

u  convection 

w  wall 

Other  symbols 

time-averaged 

oo  ambient  condition  or  infinity 

+  infinitely  small  incremental 


In  this  study,  the  reciprocating  flow  created  by  a  periodic  rever¬ 
sal  of  the  air  flow  was  investigated  in  an  effort  to  mitigate  the 
inherent  temperature  gradient  problem  of  the  conventional  battery 
system  using  the  uni-directional  coolant  flow  [1,6].  The  reciprocat¬ 
ing  flow  has  been  used  for  super-adiabatic  combustion  and  heat 
recovery  applications  [7-10].  The  reciprocating  cooling  flow  can 
be  created  by  using  flip  door  valves  and  the  unique  duct  design 
as  shown  in  Fig.  1  as  an  example  using  the  conventional  uni¬ 
directional  air  blower.  The  reversible  blower  could  be  also  used 
to  generate  the  reciprocating  flow  using  the  existing  duct  system 
for  the  uni-directional  flow  system. 

This  paper  presents  the  numerical  results  from  both  (i)  a  two- 
dimensional  computational  fluid  dynamics  (CFD)  model  created 
by  the  commercial  software,  ANSYS  FLUENT,  and  (ii)  a  lumped- 
capacitance  thermal  model  of  battery  cells  and  a  flow  network 
model.  The  battery  heat  generation  was  approximated  by  uni¬ 
form  volumetric  joule  and  reversible  (entropic)  losses  in  battery 
cells  [2,11-15].  The  accuracy  of  the  results  of  the  CFD  models  was 
validated  using  the  experimental  data  of  friction  coefficients  and 


Flow  direction  from  right  to  left 


Flow  direction  from  left  to  right 


Fig.  1.  Schematic  view  of  a  battery  system  using  a  reciprocating  air  flow  for  battery 
thermal  management.  The  flow  directions  are  (a)  from  right  to  left  side  and  (b)  from 
left  to  right  side. 
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Fig.  2.  (a)  Battery  system  in  in-line  cell  arrangement  showing  the  computational 
domain  (shaded  region)  of  a  representative  battery  submodule  with  eight  cells  con¬ 
sidered  for  the  numerical  analysis  and  (b)  the  mesh  system  of  the  battery  submodule. 


heat  transfer  coefficients  for  the  in-line  tube-bank  systems  [  1 6-2 1  ] 
which  approximate  the  battery  system  considered  for  this  study.  A 
parametric  study  using  various  operating  and  design  variables  such 
as  reciprocation  period,  cell  spacing,  coolant  flow  rate  and  different 
battery  loading  was  carried  out  and  the  results  are  discussed  in  the 
following  section. 


2.  Numerical  computation  methodology 

2.1.  Computational  domain  of  the  battery  system 

The  cylindrical  Li-ion  (LiMn204/C,  capacity:  3.6  Ah, 
D  x  L  =  42.4  mm  x  97.7  mm)  cell  was  considered  for  this  study. 
Two  methods  were  used  for  the  numerical  analysis:  (i)  two- 
dimensional  CFD  model  which  provides  precise  results  as  the 
benchmark  solutions  but  at  the  cost  of  an  intensive  computation 
(8-24  h  per  case);  (ii)  lumped-capacitance  thermal  model  for 
battery  cells  and  flow  network  model  assuming  each  battery  cell 
to  have  a  uniform  temperature  and  the  air  flow  over  the  in-line 
cells  at  constant  surface  temperatures  which  provides  reasonably 
accurate  results  on  a  light  computation  (15  min  per  case).  The 
battery  system  considered  for  this  study  consists  of  multi  rows 
(submodules)  of  a  set  of  battery  cells  as  shown  in  Figs.  1  and  2(a). 
The  shaded  region  of  the  eight-cell  submodule  as  shown  in  Fig.  2(a) 
was  used  as  the  computation  domain  for  the  numerical  analysis. 
The  symmetric  planes  were  chosen  between  the  battery  submod¬ 
ules  instead  of  the  centerline  of  the  battery  submodule.  This  is 
because  Karman  vortexes  are  created  behind  the  last  cells  at  high 
Reynolds  number  conditions  causing  an  asymmetrical  flow.  The 
computation  domain  also  includes  an  extra  flow  region  extended 
from  the  battery  area  to  avoid  developing  velocity  conditions 
affecting  the  computation  results.  The  total  length  (Zt)  of  the 
computational  domain  of  the  battery  system  in  the  streamwise 
direction  is  1060  mm  including  the  portion  of  the  battery  cells 
(lb  =  41 3.4  mm).  The  length  of  the  extension  to  the  left  side  of  the 
battery  area  is  Z*  =  323.3  mm  and  the  length  of  the  extension  to  the 
right  side  is  Z0  =  323.3  mm.  The  width  (w)  of  the  battery  submodule 
is  53  mm.  The  battery  diameter  D  is  42.4  mm.  The  dimensions  and 
specifications  of  the  battery  system  considered  for  this  study  are 
listed  in  Table  1.  Only  radial  conduction  in  the  battery  cells  was 
considered  for  this  study  because  of  the  layered-film  construction 
of  cylindrical  battery  cells  [22].  For  example,  the  conduction  ther¬ 
mal  resistance  in  the  radial  direction  of  a  spirally  wound  Li/SOCQ 
cell  would  be  twenty  times  higher  than  the  axial  direction  [23]. 


2.2.  Two-dimensional  CFD  model 

The  two-dimensional  CFD  model  was  developed  by  using  a  com¬ 
mercial  package,  ANSYS  FLUENT  12.1.4  [24]  which  is  based  on 
the  finite  volume  method  [25].  For  the  CFD  analysis,  the  buoy¬ 
ancy  effect  was  neglected;  the  Reynolds  stress  and  renormalization 
group  (Z<-£)  was  used  for  the  turbulent  model;  the  Navier-Stokes 
momentum  equation  was  solved  by  using  second-order  upstream 
difference  scheme;  the  PISO  (pressure-implicit  with  splitting  of 
operators)  algorithm  for  the  unsteady  problem  and  the  PRESTO 
(pressure  staggering  option)  for  the  swirling  flow  were  considered. 
The  reciprocating  flow  was  created  by  using  a  journal  file  available 
from  the  FLUENT  to  periodically  interchange  the  boundary  con¬ 
ditions  between  the  inlet  and  outlet,  where  a  velocity  boundary 
condition  was  used  for  the  inlet  and  a  pressure  boundary  condition 
was  used  for  the  outlet  every  reciprocation  period. 

The  shaded  region  of  the  eight-cell  submodule  as  shown  in 
Fig.  2(a)  was  used  as  the  computation  domain  for  the  CFD  analysis. 
The  mesh  systems  were  generated  by  using  commercial  meshing 
software,  GAMBIT  [26].  The  computational  domain  was  meshed 
using  structured  quadrilateral  cells  of  high  skewness  using  a  size 
function  to  assign  finer  meshes  along  the  battery  surface.  The  mesh 
quality  was  evaluated  by  checking  the  dependency  of  the  results  of 
the  CFD  analysis  such  as  pressure  drop  of  the  coolant  flow  and  maxi¬ 
mum  cell  temperature  on  various  mesh  sizes.  Finer  mesh  was  used 
for  the  flow  region,  especially  around  the  battery  cells,  whereas 
coarse  mesh  was  used  inside  the  battery  cell  volume.  Considering 
the  accuracy  of  the  results  and  computation  time,  a  mesh  system 
consisting  of  113,470  mesh  elements  for  the  baseline  system  as 
shown  in  Fig.  2(b)  was  chosen  for  the  CFD  analysis.  The  green  lines 
in  the  figure  represent  the  mesh  lines  which  were  too  fine  to  clearly 
display.  (For  interpretation  of  the  references  to  color  in  this  text,  the 
reader  is  referred  to  the  web  version  of  this  article.) 

The  computational  parameters  such  as  time  step  and  number  of 
iteration  per  time  step  used  for  the  analysis  were  chosen  to  ensure 
that  the  CFD  results  are  independent  on  the  computational  parame¬ 
ters.  The  convergence  criteria  of  the  residuals  of  the  computational 
variables  are  10-4  for  continuity  and  momentum  equations  and 
10-6  for  energy  equation.  Note  that  the  reciprocating  flow  system 
is  intrinsically  unsteady  problem  but  has  a  quasi-steady  solution 
which  is  characterized  by  a  cyclic  oscillation  every  reciprocating 
period.  The  computational  time  was  typically  about  8-24  h  to  get  a 
quasi-steady  state  solution  using  a  Window  NT  workstation  (64  bit 
OS)  with  an  Intel  Xeon  Quad-Core  processor  (3.16  GHz)  with  24  GB 
RAM. 

2.3.  Lumped-capacitance  thermal  model  and  flow  network  model 

Commercial  Li-ion  batteries  are  typically  constructed  in  a  mul¬ 
tilayer  structure  in  which  the  radial  thermal  conductivity  is  lower 
than  the  axial  one.  Nevertheless,  the  thermal  resistance  by  the 
radial  conduction  is  still  much  less  than  the  convective  thermal 
resistance  as  air  is  used  as  the  coolant  (i.e.,  Bi  =  Lchf/ks  <  0.1). 
Therefore,  a  lumped-capacitance  thermal  model  for  battery  cells 
assuming  a  uniform  temperature  in  each  cell  might  be  suffi¬ 
cient  without  compromising  accuracy  of  the  numerical  analysis 
[1-3].  The  lumped-capacitance  model,  if  properly  validated,  could 
provide  an  effective  and  simple  computational  tool  to  perform 
dynamic  thermal  analyses  using  battery  duty  cycles  and  real-time 
battery  control  on  light  and  fast  computation.  In  case  of  analyzing 
complex  battery  systems  including  duct,  battery  holder,  flow  baf¬ 
fle  and  electrical  connector/wire,  the  CFD  modeling  would  be  the 
choice  but  require  an  immense  effort  on  meshing,  and  computation 
as  compared  to  the  lumped-capacitance  modeling. 

The  battery  submodule  of  eight  cells  (shaded  region  in  Fig.  2(a)) 
was  considered  for  the  numerical  analysis  using  the  lumped- 
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Table  1 

Dimensions  of  battery  cell  and  submodule  and  the  baseline  conditions  used  for  the  numerical  analysis  and  thermophysical  properties  of  Li-ion  battery  and  air. 


Battery  technology:  Li-ion  battery  (cathode:  LiMn204,  anode:  carbon) 


Battery  cell  specification 

Air  properties  (at  300  K) 

Capacity  [Ah] 

3.6 

pf  [kgm-3] 

1.1614 

D  [mm] 

42.4 

kf  [W  m-1  K"1] 

0.0263 

L  [mm] 

97.7 

cP/  [Jkg-1  K"1] 

1007 

M  [kg] 

0.3 

M/[  kgm-1  s-1] 

1.846  x  10-5 

Aku  [mm2] 

13,014.0 

P‘7 

0.707 

Vs  [mm3] 

137,948.2 

ps  [kgirr3] 

2007.7 

ks  [Wm-1  K"1] 

32.2 

cP,s  [Jkg1  K1] 

837.4 

Baseline  conditions 

dEoddT  [mVK-1  ] 

-0.3 

Charge/discharge  rate  [C] 

7 

Battery  submodule  specification 

Tfoo  [°C] 

20 

Number  of  cell  in  flow  direction,  n 

8 

Be  D,  max 

13,300 

li  [mm] 

323.3 

Vf  [CFM] 

8.8 

l0  [mm] 

323.3 

r[s] 

1200 

lb  [mm] 

413.4 

ST  [mm] 

53 

lt  [mm] 

1060 

SL  [mm] 

53 

capacitance  cell  thermal  model.  The  thermal  circuit  diagram  of  the 
submodule  is  shown  in  Fig.  3(b)  and  (c).  The  energy  equation  for  a 
cell  (index  number  i)  can  be  given  by 

dT<  ,■ 

(pcpV)si^=SrJ-QkuJ  (1) 

Here  5r  l  is  the  heat  generation  of  a  battery  cell.  The  battery  heat 
generation  is  approximated  by  the  sum  of  (i)  Joule  (irreversible) 
and  (ii)  entropic  (reversible)  heats  and  is  given  by 

$r,i  —  Qirr  +  Qrev  (2) 

The  first  term  (Qirr)  in  the  battery  heat  generation  represents  the 
Joule  heating  based  on  the  battery  internal  electrical  resistance  Re 
and  is  given  by 

Qirr  =  I(E  ~  Eoc)  =  I2 Re  (3) 

where  the  battery  electrical  resistance  was  experimentally  deter¬ 
mined  by  averaging  over  a  range  of  the  state-of-charge  (SOC)  and 
is  expressed  in  a  function  of  the  cell  temperature  by  [2] 

Re  =  —0.0001  T3  +  0.0134 T2  - 0.5345  T  + 12.407  (4) 

where  Re  is  the  internal  electrical  resistance  [m£2  ]  of  a  Li-ion  battery 
cell.  T  is  the  cell  temperature  [°C]. 


The  second  term  (Qrev)  in  the  battery  heat  generation  represents 
the  reversible  entropic  loss  [12]  and  is  given  as 

Qrev  =  ~ IT ^  (5) 

where  the  entropy  coefficient  dE0CldT  is  a  function  of  charge  den¬ 
sity,  SOC  and  cell  temperature.  The  SOC-averaged  value  of  -0.3 
[mVK-1  ]  was  used  for  the  entropic  coefficient  dE0Cl dT  which  was 
obtained  from  the  results  of  the  experimental  measurement  of 
two  half-cells  (LiMn204/Li  and  Li/C)  [11-13].  For  the  SOC  range 
(100-70%)  considered  for  this  study  as  shown  in  Fig.  6,  the  com¬ 
position  stoichiometry  ranges  are  y  =  0.442-0.590  for  the  positive 
electrode  (LiyMn204)  and  x  =  0.676-0.511  for  the  negative  elec¬ 
trode  (Li*C6)  [14,15].  The  entropic  coefficient  of -0.3  [mVK-1  ]  was 
obtained  by  adding  the  averaged  values  of  the  dE0CldT  of  the  two 
electrodes  over  the  corresponding  stoichiometry  ranges  [11].  It  is 
worth  mentioning  that  for  hybrid  electric  vehicle  application,  the 
entropic  loss  is  often  neglected,  because  at  high  charge/discharge 
rates,  the  irreversible  Joule  heating  becomes  dominant  than  the 
reversible  entropic  heating  [14].  The  Joule  heating  quadratically 
increases  with  respect  to  electrical  current,  whereas  the  entropic 
loss  linearly  increases  with  respect  to  electrical  current  and  can  be 
either  positive  (exothermic)  or  negative  (endothermic)  depending 


Fig.  3.  (a)  Battery  submodule  with  eight  cells  used  for  the  lumped-capacitance  thermal  modeling,  (b)  Thermal  circuit  diagram  for  a  representative  battery  cell  showing  the 
battery  heat  generation  and  heat  transfer,  (c)  Thermal  network  diagram  for  the  battery  submodule  showing  the  temperature  nodes  for  battery  cells  and  fluid  stream. 
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on  the  charging  or  discharging  conditions.  As  a  result,  the  net  heat¬ 
ing  effect  of  the  reversible  heating  could  be  negligible  after  a  full 
cycle  of  charging/discharging. 

As  shown  in  Fig.  3(b)  and  (c),  the  surface  convection  heat  transfer 
rates  Qj<u  i  of  a  cell  (i  =  1,  2, . . .,  8  from  the  left  and  to  the  right  side 
as  shown  in  Fig.  3(c))  can  be  expressed  using  convection  thermal 
resistances  Ru  i  [10,21]  and  depending  on  the  flow  direction,  are 
given  by 


drop  (friction  factor)  [17,19]  and  heat  transfer  (Nusselt  number) 
[16,20,21]. 

Fig.  4  shows  the  comparison  of  the  friction  factor  results  from  the 
2D  CFD  simulation  and  the  experimental  measurement  of  the  in¬ 
line  tube-bank  systems  under  an  isothermal  condition.  The  friction 
factor  (f)  is  defined  by 


Ql<u,i  < 


'  Ts,i  ~  Tfj_  i 

Ru,i 

Ts,i  —  Tfj 
Ru,i 


for  the  coolant  flowing  from  the  left  to  the  right  side  (notated  as  solid  arrows  in  Fig.  3(b)  and  (c)) 
for  the  coolant  flowing  from  the  right  to  the  left  side  (notated  as  dotted  arrows  in  Fig.  3(b)  and  (c)) 


(6) 


where  the  convection  thermal  resistances  for  the  cells  are  given  by 
1 


Ru,i  = 


NTU  = 


(pcpV)f{  1-e-NTU) 

 \u,i  Nu Dkf 


(, PCpV)f  ( pcpv)f 


D 


(7) 


(8) 


where  the  Nusselt  number  (NuD)  proposed  by  Zukauskas  and  Ulin- 
skas  [16,21]  for  the  in-line  tube-banks  system  is  given  by 


Nud  =  < 


0-8Re°AmaxPro  :ib  (Pr/Prw)u 


^0.25 

0.51  Re°D5maxPr°-36(Pr/Prwf25 


yO-36  ( 
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103  <  keD,max  <  2  X  10b 
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0.021Re 


,0-84  Dl-0.4 


PrU4(Pr/Prw) 


2  X  10s  <  Reo,max  <  2  X  106 


R&D,  max  — 


pUmaxD 
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(10) 


Du  ("> 

Flere  U  is  the  free  space  velocity.  Umax  is  the  maximum  fluid  velocity 
in  the  transverse  spacing  (Sj  -  D )  between  two  vertically  adjacent 
cells.  The  Nusselt  number  in  Eq.  (9)  is  valid  only  when  the  number 
of  cells  in  the  streamwise  direction  exceeds  20.  For  the  eight-cell 
system  used  for  this  study,  a  correction  factor  C=0.95  [21]  is  used 
as  below 


Nup^c  =  C  Nu£ 


(12) 


f  = 


A  p 


(1/2  )npU2mi 


(15) 


where  A p  is  the  total  pressure  drop  of  the  tube-bank  system,  n  is 
the  total  number  of  the  tubes  (or  cells)  in  the  streamwise  direction. 
It  is  shown  from  Fig.  4  that  there  is  a  noticeable  transition  in  the 
friction  factor  curves  around  the  Reynolds  number  of  5  x  1 03  where 
the  flow  transition  from  laminar  to  turbulent  regimes  occurs.  It  is 
also  shown  that  the  result  from  the  2D  CFD  simulation  is  in  good 
agreement  (within  5%  variation  at  the  baseline  Reynolds  number 
ReD,max  =  13300)  with  the  experimental  results  despite  the  fact  that 
the  bounding  walls  could  introduce  a  three-dimensional  end  effect 
into  the  experimental  results. 

The  Nusselt  number  correlations  from  the  experimental  mea¬ 
surement  of  the  inline-tube-bank  systems  [16,18-21]  were  also 
compared  with  the  results  from  the  2D  CFD  simulation  as  shown 
in  Fig.  5.  Two  empirical  correlations  of  the  Nusselt  number  were 
considered  for  the  comparison  for  the  model  validation:  one  by 
Zukauskas  and  Ulinskas  [Eqs.  (9)-(l  1 )]  using  the  Reynolds  number 
based  on  the  maximum  velocity  (Umax)  in  the  least  free  space;  the 
other  by  Whitaker  [20]  using  the  Reynolds  number  based  on  the 
free  space  velocity  (D)  and  the  particle  diameter  (Dp)  as  defined 


below. 

(16) 

(17) 

The  comparison  of  the  empirical  correlation  (Eq.  (9))  and  the  results 
from  2D  CFD  modeling  is  discussed  in  Section  3. 

As  shown  in  Fig.  3(b),  the  energy  equation  for  the  fluid  stream 
flowing  over  a  cell  (i  =  1,  2, . . .,  8)  can  be  given  by 

Qu,i-i  -  Qu,i  +  Qku,i  =  o  (13) 

where 

Qu,i-1  -  Qu,i  =  (p  cpV)f(TfJ_,  -  TfJ)  =  (Mcp)f(Tf  i_-i  -  TfJ)  (14) 

The  first-order  ordinary  differential  equations  (ODEs)  for 
eight  cells  (Eq.  (1))  were  solved  using  an  initial-value  solver 
based  on  the  backward  differentiation  formulas  (BDF)  method 
[27]. 

3.  Results  and  discussion 

3. 1.  Validation  of  CFD  model  using  uni-directional  flow  system 

The  battery  cell  arrangement  used  for  this  study  can  be  approx¬ 
imated  by  the  in-line  tube-bank  system  [16-21,28]  as  shown  in 
Fig.  2.  The  2D  CFD  model  was  validated  using  the  experimental 
results  of  the  tube-bank  systems  using  a  uni-directional  flow.  The 
experimental  results  used  for  the  validation  include  the  pressure 


where  s  is  the  porosity  (or  void  volume  ratio)  of  the  portion  ( lb  as 
shown  in  Fig.  2(a))  of  the  battery  system  occupied  by  the  battery 
cells. 


Fig.  4.  Comparison  of  friction  factor  versus  Reynolds  number  between  the  CFD 
results  and  experimental  results  of  the  in-line  tube-bank  system  using  isothermal 
wall  condition  and  uni-directional  flow. 
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Fig.  5.  Comparison  of  Nusselt  number  versus  Reynolds  number  between  the  CFD  Fig.  6.  Temporal  variation  of  the  heat  generation  rates  (reversible,  irreversible  and 

results  and  experimental  results  of  the  in-line  tube-bank  system  using  isothermal  total)  for  the  battery  cell  #4  (refer  to  Fig.  3(a)),  electrical  current  and  state-of-charge 

wall  condition  and  uni-directional  flow.  (SOC)  for  the  baseline  condition. 


The  Nusselt  number  using  the  results  of  the  CFD  simulation  is 
calculated  using  a  log-mean  temperature  difference  (LMTD)  [21] 
and  is  given  by 


kt  hD 

Nud  =  — 


where 


h=  %Cp,/(r/,L  -  Tf'0) 

■AsTlmTD 


(18) 


(19) 


As  =  rntDL 


(20) 


Tlmtd  = 


(.Ts,8  -  TftL)  -  (TS't  —  Tf  o) 
M(Ts,8-Tf,L)/(TsA-TL0)) 


(21) 


Here  D  is  the  diameter  of  the  battery  cell.  L  is  the  length  of  the 
cell,  n  is  the  total  number  of  the  cells.  Ts  i  and  Tsg  are  the  temper¬ 
ature  of  the  first  and  eighth  cells,  respectively.  Tfja  and  Tf  i  are  the 
fluid  temperatures  at  the  left  and  right  sides  of  the  battery  system, 
respectively  (refer  to  Fig.  3(c)). 

As  shown  in  Fig.  5,  the  Nusselt  number  from  the  2D  CFD  analysis 
is  in  good  agreement  (within  3.5%  variation  at  the  baseline  Reynolds 
number  ReD  max  =  13,300)  with  the  experimental  results  of  the  in¬ 
line  tube-bank  system  using  the  uni-directional  flow  [16,18-21]. 
Therefore,  it  is  rational  to  conclude  from  the  comparisons  of  the 
uni-directional  flow  system  in  Figs.  4  and  5  that  the  2D  CFD  model 
is  reasonably  accurate  to  be  used  for  the  analysis  using  the  recip¬ 
rocating  flow  system. 


3.2.  Results  of  reciprocating  flow  system  using  baseline  conditions 

The  operation  conditions  and  design  variables  such  as  reciproca¬ 
tion  period,  cell  spacing,  different  battery  loading,  and  air  flow  rate 
were  considered  for  the  parametric  study  to  investigate  the  effect 
of  the  variables  on  battery  cooling  performance,  especially  battery 
cell  temperature.  The  dimensions  and  thermo-physical  properties 
of  the  Li-ion  battery  cell  and  battery  system  used  for  this  study  are 
listed  in  Table  1.  Few  important  baseline  conditions  used  for  the 
parametric  study  are  listed  below. 

•  Battery  technology:  Li-ion  battery  (cathode  material:  LiMn204, 
anode  material:  carbon;  cell  capacity  =  3.6  Ah;  diameter 
D  =  42.4  mm  and  length  L  =  97.7  mm). 

•  Battery  duty  cycle:  the  constant  stepwise  charge  and  discharge 
of  7C  (25.2  A)  with  each  half  cycle  period  of  150  s  were  used. 


The  battery  duty  cycle  was  chosen  to  maintain  the  SOC  value 
higher  than  70%  (refer  to  Fig.  6). 

•  Inlet  air  temperature:  Tft00  =  20  °C.  The  air  as  the  coolant  fluid  was 
assumed  to  be  conditioned  by  vehicle  air-conditioning  system. 

•  Reciprocation  period:  r  =  1200  s.  The  cyclic  period  (r)  of  the  recip¬ 
rocating  flow  is  defined  as  the  time  for  a  reciprocating  flow  to 
recover  its  initial  flow  direction.  Note  that  the  uni-directional 
flow  with  no  change  in  flow  direction  is  a  special  case  of  the 
reciprocating  flow  having  an  infinite  reciprocation  period  (r  =  oo). 

•  Battery  cell  spacing:  ST/D  =  1.25  was  used  for  the  transverse  direc¬ 
tion  and  SlID  =  1.25  was  used  for  the  longitudinal  (streamwise) 
direction.  Note  that  S/D  =  1.0  means  no  gap  between  the  cells 
(refer  to  Fig.  2(a)). 

•  Air  flow  rate:  the  volumetric  flow  rate  Vf  =  8.8  CFM  per  a  battery 
submodule  of  eight  cells  was  used  which  is  equivalent  to  the  fluid 
velocity  of  1  m  s-1  or  Reo.max  =  13,300. 

Fig.  6  shows  the  profiles  of  the  electrical  current,  SOC,  and 
battery  heat  generation  of  the  fourth  cell  (#4)  chosen  as  a  repre¬ 
sentative  cell  of  the  battery  system  (refer  to  Fig.  3).  The  battery 
heat  generation  considered  for  this  study  consists  of  irreversible 
(Joule)  and  reversible  (entropic)  parts  as  defined  in  Eq.  (2).  The  irre¬ 
versible  part  is  always  exothermic  and  dominant  for  higher  current 
conditions.  Whereas,  the  irreversible  part  can  be  either  exothermic 
or  endothermic,  thus  it  could  be  canceled  out  after  a  full  cycle  of 
charging/discharging. 

Fig.  7  shows  the  results  of  the  temperature  and  flow  veloc¬ 
ity  from  the  CFD  analysis  using  the  baseline  conditions.  Fig.  7(a) 
shows  the  temperature  contour  of  the  battery  submodule  at  various 
times  within  the  reciprocating  period  of  r  =  1 200  s.  The  temperature 
results  were  obtained  from  the  6th  cycle  (t  =  6000-7200  s)  when  a 
quasi  steady-state  was  achieved.  The  quasi  steady-state  is  defined 
as  the  state  where  the  oscillation  amplitudes  of  the  cell  tempera¬ 
tures  converge  to  constant  values.  The  quasi  steady-state  condition 
is  more  discussed  in  the  next  section.  The  direction  of  the  recipro¬ 
cating  flow  at  the  beginning  of  the  6th  cycle  tc(  =  t-  4r)  =  0+  is  from 
the  right  to  the  left  side.  Shortly  after  that  moment  (within  few 
milliseconds),  the  flow  direction  is  reversed,  i.e.,  from  the  left  to 
the  right  side  and  quickly  becomes  fully  developed.  At  the  end  of 
the  half  cycle  tc  =  r/2,  the  flow  direction  is  reversed  again.  At  the 
end  of  the  full  cycle  tc  =  r,  the  temperature  profile  is  recovered  to 
the  initial  condition  at  tc  =  0+.  It  is  also  observed  from  Fig.  7(a)  that 
the  temperatures  of  a  pair  of  the  outer  cells  located  downstream  of 
the  flow  (i.e.,  #1  and  2  at  tc  =  0+  and  #7  and  8  at  tc  =  r/2)  are  close 
each  other.  From  a  close  examination  of  the  velocity  results,  it  is 
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Fig.  7.  Temporal  progression  of  (a)  temperature  contours,  (b)  battery  cell  temperatures,  and  (c)  velocity  vector  contours  of  the  battery  submodule  under  the  baseline  condition 
(r  =  1200s). 


observed  that  it  takes  only  few  milliseconds  for  the  old  flow  pat¬ 
tern  to  fade  away  and  for  the  new  flow  pattern  to  fully  establish 
when  the  flow  direction  is  reversed. 

Fig.  7(b)  shows  the  temporal  variation  of  the  area-averaged  tem¬ 
perature  of  each  cell  at  various  times.  It  is  observed  from  the  figure 
that  the  temperatures  pivot  at  the  center  (between  #4  and  5  cells) 
of  the  battery  system.  As  a  result,  the  temperatures  of  the  inner 
cells  vary  less,  while  the  outer  cell  temperatures  oscillate  more. 


The  figure  also  shows  that  the  maximum  (both  time-averaged  and 
instantaneous)  temperatures  of  the  reciprocating  flow  are  lower 
than  the  unidirectional  case  (r  =  oo).  Especially,  the  maximum  dif¬ 
ference  of  the  time-averaged  temperatures  (notated  in  red  dotted 
line)  is  1.5  °C  which  is  much  lower  than  the  uni-directional  case 
of  about  5.5  °C.  It  is  worth  noting  that  the  instantaneous  max¬ 
imum  temperature  occurs  at  the  second  (#2)  and  seventh  (#7) 
cells,  while  the  time-averaged  maximum  temperature  occurs  at 
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Fig.  8.  Temporal  variations  of  the  battery  cell  temperatures  for  various  reciprocation  periods  and  the  total  battery  heat  generation  of  the  battery  cell  #4  (refer  to  Fig.  3(a)): 
(a)  r  =  oo,  (b)  r  =  1200  s,  and  (c)  r  =  120  s.  (d)  A  magnified  view  of  the  quasi  steady-state  temperatures  for  the  reciprocation  periods  of  r  =  120  s. 


the  cells  (#4  and  5)  near  or  at  the  middle  of  the  battery  system. 
The  time-averaged  maximum  temperature  is  the  maximum  of  the 
time-averaged  temperature  of  the  cells  over  a  reciprocating  period. 
The  definition  of  the  time-averaging  of  the  cell  temperatures  is 
discussed  in  the  next  section. 

Fig.  7(c)  shows  the  velocity  vector  field  at  various  times  within 
the  reciprocating  period  of  r  =  1200  s.  In  the  figure,  it  is  observed 
that  the  air  flow  accelerates  and  decelerates  while  passing  through 
the  least  gap  (the  transverse  spacing  ST  -  D )  between  the  cells  and 
the  flow  separations  are  created  after  each  cell,  except  that  rather 
large  recirculation  is  created  behind  the  last  cell  with  respect  to 
the  flow  direction.  It  is  also  observed  that  the  stream  through  the 
transverse  spacing  is  prevailing  than  one  through  the  longitudinal 
spacing  (S l  -D).  This  flow  channeling  effect  is  severer  with  nar¬ 
rower  longitudinal  spacing  and  is  very  detrimental  for  the  battery 
cooling  because  of  ineffective  use  of  the  battery  surface  for  forced 
convective  heat  transfer,  which  could  be  mitigated  by  using  a  larger 
longitudinal  spacing  and  flow  deflectors. 


3.3.  Effect  of  reciprocation  period 

The  temporal  variations  of  the  battery  cell  temperatures  and 
heat  transfer  rates  for  various  reciprocation  periods  are  shown 
in  Figs.  8  and  9.  The  temperature  and  heat  transfer  results  were 
obtained  from  the  2D  CFD  model  for  three  different  reciprocating 
periods:  r  =  oo,  1 200  and  1 20  s. 


Figs.  8  and  9(a)  show  the  results  of  the  “uni-directional  flow 
case”  which  is  a  special  case  of  the  reciprocating  flow  with  r  =  oo 
as  the  coolant  constantly  flows  from  the  left  to  the  right  side.  In 
Figs.  8  and  9(a),  a  periodic  fluctuation  of  the  temperature  and  heat 
transfer  rate  is  shown  due  to  the  alternating  sign  of  the  entropic 
heat  (Qrev)  depending  on  charging/discharging  of  the  battery  cycle 
with  a  period  of  150  s  as  shown  in  Fig.  6.  Flowever,  the  temper¬ 
atures  and  heat  transfer  rates  quickly  converge  and  approach  to 
quasi  steady-states  around  5000  s  as  shown  in  Figs.  8  and  9(a). 
The  time-averaged  temperatures  were  calculated  after  reaching  the 
quasi  steady-state  over  a  period  (300  s)  of  the  battery  duty  cycle. 
Note  that  since  the  heat  generation  is  not  constant  due  to  the  sign- 
alternating  entropic  heat,  the  maximum  temperature  and  average 
maximum  temperature  are  different  for  the  uni-directional  flow 
case  as  shown  in  Fig.  10. 

It  is  also  observed  from  Fig.  8(a)  that  the  lowest  temperature 
occurs  at  the  first  cell  (#1)  and  the  highest  temperature  occurs  at 
the  seventh  cell  (#7)  and  the  maximum  cell  temperature  differ¬ 
ence  ( ATs.max)  is  5.5  °C  for  the  uni-directional  flow  case.  Although 
the  highest  cell  temperature  is  expected  to  occur  at  the  last  cell  (#8) 
near  to  the  flow  exit  due  to  the  highest  air  temperature,  the  highest 
cell  temperature  occurs  at  the  seventh  cell  (#7).  This  is  because  the 
eighth  cell  is  exposed  to  more  convective  environment  due  to  big¬ 
ger  recirculation  formed  behind  the  cell  as  shown  in  Fig.  7(c)  than 
the  inner  cells.  It  is  worth  mentioning  that  the  lumped  thermal 
model  using  the  empirical  correlations  (Eqs.  (l)-(8))  always  pre¬ 
dicts  the  highest  temperature  in  the  last  cell  (#8)  which  attributes 
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Fig.  9.  Temporal  variations  of  the  heat  transfer  rates  for  various  reciprocation  periods  and  the  total  battery  heat  generation  of  the  battery  cell  #4  (refer  to  Fig.  3(a)):  (a)  r  = 
(b)  r  =  1200  s,  and  (c)  r  =  120  s.  (d)  A  magnified  view  of  the  quasi  steady-state  heat  transfer  rates  for  the  reciprocation  periods  of  r  =  120  s. 


to  the  fact  that  the  empirical  correlations  predicts  the  averaged  heat 
transfer  for  the  cells  in  the  system. 

For  the  “reciprocating  flow  cases”  as  shown  in 
Figs.  8(b)-(d)  and  9(b)-(d),  the  coolant  flow  direction  is  reversed 
at  every  half  cycle  (tc  =  r/2)  of  each  reciprocating  period.  Fig.  8(b) 
and  (c)  shows  the  cell  temperatures  and  the  total  battery  heat 
generation  of  the  fourth  cell  (#4)  for  two  different  reciprocating 


Fig.  10.  Effect  of  the  reciprocation  period  on  the  cell  temperatures  (instantaneous 
maximum:  time-averaged  maximum  and  minimum  and  differential). 


periods  of  r  =  1 200  and  1 20  s.  Fig.  9(b)  and  (c)  shows  the  convective 
heat  transfer  rates  of  the  cells  and  the  total  battery  heat  generation 
of  the  fourth  cell  (#4)  for  the  corresponding  reciprocating  periods. 
Figs.  8(d)  and  9(d)  show  the  magnified  view  of  the  inset  boxes 
shown  in  Figs.  8(c)  and  9(c). 

Fig.  8(b)  shows  that  the  temperatures  of  the  outer  cells  (#1  and 
#8)  near  to  the  flow  inlet  and  outlet  have  a  bigger  swing  than  the 
inner  cells  due  to  the  larger  change  in  the  ambient  temperature 
and  strong  convective  heat  transfer  by  the  jet  impingement  of  the 
approaching  flow  which  is  also  shown  in  Fig.  9(b).  The  symmetric 
oscillatory  behavior  of  the  cell  temperatures  as  shown  in  Fig.  7(b) 
happens  for  all  reciprocating  periods.  It  is  observed  that  the  larger 
the  reciprocating  period  is,  the  bigger  the  swings  in  the  cell  temper¬ 
atures  and  convective  heat  transfer  rates  are.  This  is  because  of  the 
longer  time  for  the  cell  temperature  to  change  resulting  in  the  larger 
temperature  difference  between  the  cells  and  air  temperature.  Note 
that  regardless  of  the  reciprocating  period,  the  amplitudes  of  the 
temperature  swing  tends  to  converge  to  quasi  steady-states  around 
5000  s,  similarly  found  in  the  uni-directional  case. 

For  relatively  shorter  reciprocating  periods  (r  =  120  s)  as  shown 
in  Fig.  8(c)  and  (d),  the  temperatures  of  the  inner  cells  are  always 
higher  than  the  outer  cells.  This  is  due  to  the  heat  accumulation 
effect  by  the  fast  reversal  of  the  reciprocating  flow.  From  the  prac¬ 
tical  viewpoint  of  battery  cooling  system  operation,  smaller  period 
would  be  not  desirable  for  actual  applications. 

Fig.  10  shows  the  variation  of  the  time-averaged  temperatures 
with  the  reciprocation  period,  which  were  calculated  over  a  full 
reciprocating  period  after  the  quasi  steady-state  is  achieved.  The 
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Fig.  11.  Effect  of  the  battery  spacing  on  the  maximum  cell  temperature  and  pressure 
drop. 


maximum  (  Ts,max)  and  minimum  (Ts?min)  time-averaged  tempera¬ 
tures  are  defined  by 


Ts,  max  — 


Ts,  max  dt 
X 

Ts,  min  dt 
X 


(22) 

(23) 


where  Ts,max  could  indicate  the  level  of  the  thermal  stress  (or 
health)  of  a  battery  system.  The  difference  between  the  maximum 
and  minimum  of  the  time-averaged  cell  temperatures  is  given  by 


A  Ts. 


Jo  (?s,  max  —  Ts  m jn)  dt 


=  TS. 


-  T- 


(24) 


ATs^max  is  related  to  the  level  of  the  imbalance  of  the  thermal  stress 
of  a  battery  system. 

It  is  shown  in  Fig.  10  that  as  the  reciprocating  period  decreases 
7s, max  decreases  and  Ts  m in  increases,  thus  A Ts,max  decreases 
(i.e.,  more  uniform  temperature).  Note  that  the  instantaneous 
maximum  cell  temperature  TS)max  is  always  higher  than  the  time- 
averaged  one  Ts, max-  It  is  true  for  the  uni-directional  case  (r  =  oo) 
because  of  the  temperature  swing  due  to  the  sign-alternating 
entropic  heat.  It  is  also  shown  in  Fig.  10  that  the  maximum  of  the 
time-averaged  temperatures  Ts?max  for  the  reciprocating  flow  case 
of  r  =  120  s  is  28.1  °C  while  rs?max  =  29.6  °C  for  the  uni-directional 
flow  (r  =  oo).  A Ts,max  for  the  reciprocating  flow  of  r  =  120s  is 
significantly  reduced  to  1.5  °C  from  A Ts,max  =  5.5  °C  for  the  uni¬ 
directional  flow  case. 


3.4.  Effect  of  battery  cell  spacing 

The  cell  spacing  is  defined  in  both  horizontal  (longitudinal  spac¬ 
ing  Sl)  and  vertical  directions  (transverse  spacing  St)  as  shown  in 
Fig.  2(a).  Fig.  11  shows  the  effect  of  the  cell  spacing  on  the  max¬ 
imum  cell  temperature  of  the  battery  system  under  the  baseline 
conditions.  The  maximum  cell  temperature  is  the  instantaneous 
value  of  the  maximum  cell  temperature  in  the  battery  system  at 
any  given  time.  As  shown  in  Fig.  7(b),  the  maximum  cell  tempera¬ 
ture  always  occurs  at  the  second  (#2)  and  seventh  (#7)  cells  at  the 
time  of  the  flow  reversal  (tc  =  0  or  r/2)  because  of  the  increased  heat 
transfer  of  the  most  outer  cells  (#1  and  #8)  due  to  the  large  flow 
recirculation  behind  them. 

It  is  observed  from  Fig.  1 1  that  as  the  transverse  spacing  (Sr) 
decreases,  the  cell  temperature  are  decreased  because  the  flow 
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Fig.  12.  Effect  of  the  battery  loading  and  air  flow  rate  on  the  maximum  cell  temper¬ 
ature. 

velocity  is  increased  due  to  the  decreased  flow  area  resulting  in 
increased  convective  heat  transfer  and  pressure  drop.  It  is  also 
observed  from  the  figure  that  longitudinal  spacing  (Sl)  gets  larger, 
the  maximum  cell  temperature  is  decreased,  this  is  because  the 
flow  creates  bigger  eddies  in  the  longitudinal  spacing  of  the  cells 
rather  than  just  passing  over  the  cells  (flow  channeling)  as  shown 
in  Fig.  7(c).  As  a  result,  the  pressure  drop  of  the  system  is  increased. 
Therefore,  it  is  summarized  that  the  cell  temperature  is  decreased 
with  larger  longitudinal  and  smaller  transverse  spacing  and  vice 
versa  for  the  pressure  drop.  This  competing  result  suggests  the 
existence  of  the  optimum  spacing  for  the  lowest  maximum  cell 
temperature  for  a  fixed  pumping  power  consumption  of  the  air 
flow.  The  flow  channeling  as  previously  discussed  in  Fig.  7  could 
be  mitigated  by  using  a  flow  baffle  which  is  commonly  used  to 
redirect  the  flow  in  commercial  battery  systems.  Note  that  smaller 
spacing  would  be  preferred  because  of  the  packaging  constraints 
of  the  battery  systems  in  vehicles. 

3.5.  Effect  of  different  battery  loading 

Fig.  12  shows  the  effect  of  the  battery  charge/discharge  rate 
on  the  instantaneous  maximum  cell  temperature  under  the  base¬ 
line  conditions.  For  this  calculation,  the  battery  charge/discharge 
in  the  stepwise  profile  as  shown  in  Fig.  6  was  used  by  varying  only 
the  magnitude  of  the  electrical  current.  The  battery  heat  genera¬ 
tion  was  calculated  by  the  sum  of  Joule  (irreversible)  and  entropic 
(reversible)  heats  as  given  in  Eq.  (5).  The  reversible  entropic  heat¬ 
ing  has  a  linear  dependence  with  the  electrical  current,  thus 
the  net  heating  effect  is  zero  during  a  full  cycle  of  the  charg¬ 
ing/discharging.  The  irreversible  Joule  heating  is,  however,  always 
positive  because  of  the  quadratic  dependence  with  the  current  and 
becomes  more  dominant  in  the  battery  heat  generation  at  higher 
charging/discharging  rates.  As  shown  in  Fig.  12,  the  maximum  cell 
temperature  increases  quadratically  with  the  charge/discharge  rate 
because  of  the  non-zero  Joule  heating  effect  during  the  current 
cycling. 

Fig.  12  also  shows  the  effect  of  the  air  flow  rate  on  the  instan¬ 
taneous  maximum  cell  temperature.  The  instantaneous  maximum 
temperature  always  occurs  at  the  second  cell  (#2)  or  at  the  seventh 
cell  (#7)  of  the  battery  system  depending  on  the  flow  direction  as 
previously  discussed  with  Fig.  7.  It  is  observed  from  Fig.  12  that 
the  lower  air  flow  rate  (ReD  max  =  13,300)  is  insufficient  to  limit  the 
cell  temperature  below  40  °C  at  the  high  loading  condition  (>1 0  C). 
Fig.  12  also  shows  the  comparison  of  the  instantaneous  maxi- 


R.  Mahamud,  C.  Park  /  Journal  of  Power  Sources  196(2011)  5685-5696 


5695 


Fig.  13.  Effect  of  the  air  flow  rate  on  the  Nusselt  number  and  pressure  drop. 

mum  temperatures  from  the  CFD  model  and  a  lumped-capacitance 
model.  The  results  from  the  CFD  and  lumped  thermal  models 
were  well  matched,  especially  for  the  lower  Reynolds  number 
(ReD,max  =  13,300). 

3. 6.  Effect  of  air  flow  ra  te 

As  the  air  flow  rate  is  increased,  the  convective  heat  transfer  on 
the  cells  is  expected  to  be  increased,  thus  the  cell  temperatures  are 
decreased  but  at  the  expense  of  higher  pressure  drop  (or  pump¬ 
ing  power).  Figs.  13  and  14  show  the  effect  of  the  air  flow  rate  on 
the  heat  transfer  and  instantaneous  maximum  cell  temperature, 
respectively.  Because  of  the  transient  nature  of  the  cell  temper¬ 
atures  due  to  the  reciprocating  flow,  the  log-mean  temperature 
difference  using  the  time-averaged  cell  temperatures  over  a  full 
cycle  of  the  reciprocating  flow  was  used  for  the  calculation  of  the 
Nusselt  number  in  Fig.  13  (refer  to  Eq.  (22)). 

Fig.  13  shows  that  the  Nusselt  number  increases  more  rapidly 
until  ReD,max  =  13,300  where  the  transition  between  the  laminar 
and  turbulent  regimes  occurs  (also  refer  to  Fig.  4).  In  Fig.  13,  the  Nus¬ 
selt  number  of  the  CFD  model  for  the  reciprocating  flow  cases  was 
also  compared  with  the  empirical  correlation  [Eq.  (9)]  developed 
for  the  uni-directional  flow  system  [16,19,21].  As  shown  in  Fig.  13, 
for  the  Reynolds  numbers  below  ReD,max  =  26,600,  the  Nusselt  num¬ 
ber  of  the  CFD  model  is  greater  than  the  empirical  correlation  and 


Fig.  14.  Effect  of  the  air  flow  rate  on  the  maximum  cell  temperature. 


vice  versa  for  the  Reynolds  numbers  above  ReD,max  =  26,600.  Note 
that  the  similar  result  is  also  shown  in  Fig.  5.  Fig.  13  also  shows  that 
the  pressure  drop  rapidly  increases  quadratically  with  the  Reynolds 
number. 

Fig.  14  shows  the  change  of  the  instantaneous  maximum  cell 
temperature  as  the  Reynolds  number  is  increased.  As  shown  in 
the  figure,  the  maximum  cell  temperature  is  quickly  decreased 
until  ReD,max  =  13,300  where  the  slope  change  of  the  Nusselt  num¬ 
ber  happens  due  to  flow  transition  as  shown  in  Fig.  13.  After  the 
Reynolds  number,  the  cell  temperature  rather  slowly  decreases  as 
the  air  flow  rate  increases.  The  optimum  conditions  for  higher  heat 
transfer  and  lower  pressure  drop  (or  pumping  power)  could  exist. 
The  figure  also  shows  the  comparison  of  the  temperatures  of  the 
2D  CFD  and  lumped  thermal  models.  The  maximum  cell  temper¬ 
ature  from  the  lumped  thermal  model  is  always  lower  than  the 
CFD  model.  This  attributes  to  the  fact  that  the  empirical  correlation 
used  for  the  lumped  thermal  model  will  predict  the  averaged  heat 
transfer  for  the  entire  cells.  It  is  concluded  from  Figs.  13  and  14 
that  the  heat  transfer  correlations  for  the  uni-directional  flow  sys¬ 
tems  is  reasonably  accurate  enough  to  be  used  for  the  reciprocating 
flow  systems.  If  a  proper  heat  transfer  correlation  is  available, 
the  lumped-capacitance  thermal  model  would  be  very  suitable  for 
real-time  battery  thermal  analyses  requiring  time-efficient  and  fast 
computation. 

4.  Conclusion 

A  reciprocating  air  flow  was  used  for  battery  thermal  man¬ 
agement  of  electrical-drive  vehicles  to  improve  temperature 
uniformity,  and  reduce  maximum  cell  temperature.  The  recipro¬ 
cating  flow  for  the  battery  thermal  management  of  cylindrical 
Li-ion  (LiMn204/C)  cells  was  numerically  analyzed  using  (i)  a  two- 
dimensional  CFD  (computational  fluid  dynamics)  model  and  (ii) 
a  lumped-capacitance  thermal  model  for  battery  cells  and  a  flow 
network  model.  The  results  of  the  CFD  model  was  validated  with 
available  experimental  results  (heat  transfer  and  pressure  drop)  for 
in-line  tube-bank  systems  which  approximates  the  in-line  arrange¬ 
ment  of  the  battery  cells  considered  for  this  study.  The  numerical 
results  showed  that  the  shorter  the  reciprocating  period  is,  the 
lower  the  cell  temperature  difference  and  the  maximum  cell  tem¬ 
perature  (of  both  instantaneous  and  time-averaged)  of  the  system. 
It  was  shown  from  the  numerical  results  that  the  reciprocating 
flow  using  a  reciprocation  period  r  =  1 20  s  reduces  the  cell  temper¬ 
ature  difference  by  about  4.0  °C  (72%  reduction)  and  the  maximum 
cell  temperature  by  1.5  °C  as  compared  with  the  uni-directional 
case  (r  =  oo).  The  improvement  by  the  reciprocating  flow  mainly 
attributes  to  the  heat  redistribution  and  disturbance  of  the  bound¬ 
ary  layers  formed  on  the  cells  due  to  the  periodic  flow  reversal. 
The  results  from  the  lumped-capacitance  thermal  model  and  flow 
network  model  were  in  good  agreement  with  the  2D  CFD  results. 
Although  the  multi-dimensional  CFD  modeling  is  essential  for 
the  analysis  of  complex  battery  systems,  time-efficient  lumped- 
capacitance  models  using  appropriate  heat  transfer  correlations 
would  be  very  suitable  for  the  intensive  computations  such  as  bat¬ 
tery  lifetime  prediction  based  on  battery  duty  cycles  and  real-time 
battery  temperature  forecast  and  control. 
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